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Q_i' Abstract 

(D : 

. We study the solutions of the Einstein equations in {d + 2)-dimensions, describing 

parallel p-branes (p = d — 1) in a space with two transverse dimensions of positive gaussian 
curvature. These solutions generalize the solutions of Deser and Jackiw of point particle 
sources in (2 + l)-dimensional gravity with cosmological constant. Determination of the 
metric is reduced to finding the roots of a simple algebraic equation. These roots also 
determine the nontrivial "warp factors" of the metric at the positions of the branes. We 
discuss the possible role of these solutions and the importance of "warp factors" in the 
context of the large extra dimensions scenario. 
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1 Introduction and summary. 



We study solutions of the Einstein equations of d + 2 dimensional gravity with a cosmological 
constant. The solutions represent {d — l)-branes embedded in a space with two compact trans- 
verse dimensions with a positive curvature. They generalize the static solutions of Deser and 
Jackiw describing point particles in 2 + 1 dimensional gravity with cosmological constant, 
to the case of extended objects. 

The solutions we find obey Einstein's equations, as well as the equation of motion of the 
branes (the geodesic equation). It is remarkable that the metric, in appropriate coordinates, 
can be found without explicitly solving the Einstein equations. Analyzing the properties of the 
solutions is reduced to finding the roots of a simple algebraic equation. 

From the point of view of applications to the recently proposed [0 , large extra dimensions 
scenario, the most important property of our solutions is that they are characterized by non- 
trivial "warp factors" that affect the values of the parameters of the world- volume theories [Q. 
In particular, the usual relation between the four dimensional Planck scale, the 6 dimensional 
"fundamental" scale of gravity, and the volume of the compactified space is affected by the 
presence of the warp factors. This difference can be important if the warp factors significantly 
deviate from unity (which, in the present framework can only be achieved by fine tuning various 
parameters) (for related papers with one transverse dimension, see [§]). 

While our recipe for finding solutions is rather general, motivated by the large extra dimen- 
sions proposal, we discuss in some detail the solution that describes two 3-branes embedded in 
six dimensional space time with negative^ (de Sitter) cosmological constant. The two branes 
are located on the north and south poles of a "sphere", with a "wedge" cut out, due to the 
deficit angle characteristic of pointlike sources in 2 dimensions. An important property of the 
solution is that the metric on the branes is not fiat, but rather de Sitter. While the cosmo- 
logical constant and the strength of gravity on the "visible" brane can be tuned to satisfy the 
experimental constraints, within our ansatz there is no solution with fiat world volume metric. 

This paper is structured as follows. In Section 2, we present the ansatz and derive the 
equations of motion. We show that they reduce to a simple differential equation in one real 
variable, similar to the case considered by Deser and Jackiw In Section 3, we give a general 
analysis of the solutions. We show that the form of the metric (in appropriate coordinates) 
as well as the interpretation of its singularities can be obtained without explicitly solving the 
differential equation. All information on the solution can be deduced by finding the appropriate 
roots of an algebraic equation. In particular, these roots determine the values of the warp 
factors on the branes. We derive expressions for the strength of the gravitational coupling, the 
cosmological constant on the various branes, and the size of the transverse space. Finally, in 
Section 4, we discuss an explicit example, with = 4, which is of interest to the large extra 
dimension scenario. We show that the parameters allow fine tuning of the Planck scale and 
cosmological constant to values that do not contradict the observed ones. 



^^We note that our conventions for the signs of the curvature are those of Q; in addition A < in (|l| 
corresponds to de Sitter space. 
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2 The ansatz and equations of motion. 

The action we consider is that of gravity in D — d + 2 dimensions, with the usual Einstein 
action with cosmological term included: 

'gravity = -M^'^ J d^'y^ {R - 2A) . (1) 

The indices M,N — 1, D, while /x, i/ = 0, d ~ D — 2 and i,j = 1, 2 span the rest of the 
space; the metric has signature (—,+,...,+). The vacuum Einstein equations follow from the 
variation of the gravity action: 

3S,ra.uy = -M^'^ J d^y^g {-R^^ + \g^^R - g^^^ K) 6gMN = . (2) 

The "matter" is in the form of d — 1 dimensional branes, whose action is proportional to the 
area of the world surface they sweep in the d-\-2 dimensional spacetime: 



^matter 



Efaj d^'y I d'a S^{y-Xa{a)) ^5, (3) 



where X^{a) is the embedding of the world surface in space time (a runs over the various 
branes) and = X^ QMNiu) the induced metric and g = — det g^j^i, (hereafter we will 
omit the sum over the various sources, it will be implicit in all our formulae). The matter 
energy momentum tensor is defined by: 

SS^atter = \ j d"" y ^ T'''' {y)5gMN{y) 

= -\f j d^'y J d^a 5^iy - X{a)) g'/' X^ W(y) , (4) 

hence 

^T^N^y^ = -fjd'^a 6^{y - X{a)) f/' X^ X^ ^ • (5) 

The equation of motion of the branes (obtained by varying the action with respect to the brane 
coordinates X{a)) is the "geodesic" equation: 

where = g'^^'^d^g^^'^g^'^d^ is the covariant D'Alembertian in the induced metric. 

Our ansatz for the metric includes a nontrivial warp factor e^^, depending only on the 
coordinates transverse to the branes (denoted by x\ i — 1,2): 

ds" = e'^^"') [g^,{x'')dx^'dx'' + e^'^^^^'^dx'dx^ . (7) 

We assume that the d-dimensional metric g^ui^^) obeys the equation 

R^u - ]j^9nuR + otg^,u = , (8) 
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so that a becomes a parameter of our solution. 

The idea for finding a solution with delta-function sources is to first solve the source-free 
equations of motion. Then, admitting delta function like curvature singularities of the metric 
will allow us to include the effect of the sources. Following this strategy, in the remainder of 
this section we consider the source free equations. 

The Einstein equation following from the variation of the action with respect to the /i, v 
components of the metric is: 

,,,,,, , + ^^A,.A,^e--A-.-„ (9) 

This equation, away from the singularities, is a consequence of the other equations of motion. 
Static sources will lead to delta-function singularities on the r.h.s. of equation (^. These will 
appear, as we will show in Section 3, because of singularites of K^a. 

The (trace of) the variation of the action w.r.t. the j components of the metric gives the 
following equation: 

dd\ + ddXd\ = ^e^'+'^'A - Jl^e^^a . (10) 

We have introduced complex coordinates, z = {x^ + ix'^)/2,z = {x^ — zx^)/2, and defined 
d = dxj^ — idx2, 9 = dxj^ + idx2 such that dz = l,dz = 0. On the other hand, the traceless 
part of the equations of motion that follow from varying the action with respect to the z,j 
components of the metric are nothing but the Cauchy-Riemann conditions for the function 

V = e-2^-^aA . (11) 

The function V , therefore is holomorphic, dV = 0. This fact will play a central role in our 
ability to find explicit solutions.^ 

We can now further simplify our equations, following Multiplying (|IU]) by V given by 
( P^D and using BV = 0, we obtain after integrating over z, 

where e{z) is an integration constant. Upon introducing the new complex variable 

(13) 



V{w) 



T^^Before continuing with the analysis of solutions with nontrivial warp factors, we note that the source-free 
equations of motion admit a simple solution with K = — log(l + 4|zp/p^), and A = const., provided the 

values of A and a are related, so that the r.h.s. of ( [To| ) vanishes. These solutions are of the form dSd-2 x S'^, 
with the radius of S*^, p, and the cosmological constant of the dSd-2, cx, related by (p^: a — —2p^^{d — 2), 
A — ~2p^^d (for A = 0). These vacuum solutions of the d + 2 dimensional Einstein equations with cosmological 
constant generalize the Nariai solutions |7| to d + 2 dimensions. We thank N. Kaloper for pointing out to us 
the existence of ref . |^] . 
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as well as denoting = N, we can rewrite the equation as: 

^iV^ = J^iV^+i _ ?1 aN"-' + e . (14) 

d + 1 {d-2){d-l) ^ ' 

Now upon taking the complex conjugate equation, we see that the reality condition for 
implies that is a function of the real part of only and that e must be a real constant. 
Finally, introduce the variable 

t = TTTTTl (^ + r) = TTTT^l / TTT^ + ^.c. (15) 



{d+l)d ^ ' {d+l)d \J V{w 
and write the equation in terms of the single real variable t (p!5|): 

f{N) dN ^ ^ dN =dt, (16) 

where f{N) is defined by the first equality, and the polynomial P{N) is: 

P(N) = N'^+^ - aN'^-'^ + b , (17) 

with the coefficients a, b defined as follows: 

a d{d+l) e d+1 

^ = A id-2)id-i) ' ^ = A-r- 

This completes our discussion of the equations of motion. We saw that the equation of 
motion for the ansatz (|^) reduced to a simple differential equation in one real variable (|l^); 
a general solution can be expressed in terms of the general solution of this equation and a 
holomorphic function V (pAl). We also note that this provides also a solution to eq. — it is 
easy to see that (|TU|) together with the holomorphicity of V implies (^, away from possible 
singularities (which, as we will see in the next section, appear in K^a). 

In the following section we present a general analysis of the solutions of ([T^). 



3 General analysis of the solutions. 

In this section, we discuss the properties of the general solution of the equation (|16|). To begin, 
we need to specify the range of the variable t in ([T6|). In the cases of interest to us t will span 
the entire real axis — consider e.g. the case V{z) = z/c, so that t ~ log|2;|. In order to solve 



(16) we need to find a function F{N) such that F'{N) = f{N). Moreover, we should be able 
to invert F{N) for all values of t— then N{t) = F~'^{t) will satisfy (|l|). 

Even though the equation is rather complicated and the explicit form of F{N) is generally 
unknown, we will see that the quantities of interest of us do not require an explicit knowledge 
of the solution. In fact, the metric can, in appropriate coordinates, be expressed in terms of 
the polynomial P{N), while the ranges of the variables are determined by appropriately chosen 
roots of the polynomial P{N). 
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The basic requirement is that we find a function F{N) which monotonically changes from 
— oo to cxD over some finite range of values of N. For that to be the case, the function f{N) 
has to have the properties (recall that F' = f) that between two values of N a) it does not 
change sign and b) its modulus approaches infinity at the two boundary values of A^. Consider 
two "nearest neighbor" roots of the polynomial P{N), Ni and N2. Assume that both A^^i and 
N2 have the same sign.0 Now, at the roots of P{N) f{N) blows up; moreover, if Ni and 
N2 have the same sign, f{N) does not change sign as varies between Ni and A'2 (we are 
mostly interested in the case of even d). Then, between these two roots F{N) is monotonic 
and approaches 00 at one of the roots, say Ni, and —00 at the other root, A^2- Hence, we can 
invert the equation ([T6| ) and find the function N{t) = F^^(t) for all values of t on the real 
axis; moreover, the function N{t) approaches finite values — the roots of P{N), Ni and N2 — as 
\t\ — > 00. 

Thus the problem of finding a solution of (|16D is reduced to finding the condition that there 
are two real "nearest-neighbor" roots of P{N) of the same sign. It is easily seen, e.g. for d = 4, 
that this requires that a > 0, i.e. a and A have the same sign. 

We assume now that the parameters a, b (p!8| ) are such that two such roots, A^i and N2 are 
found (assume for simplicity that they are both positive). Let us now analyze the properties 
of the solution. The metric has the form (|^; having found the function N{t), we can now 
determine the N'^e^^ factor in the metric by using (|TT|) and the equation (|16]): 



1 

V 



dN 



2A 1 P{N) 



(19) 



Now note that using 



dt 



2A 



de = if3 



d{d+l) 
dz 



dz dz 
+ 



V{z) V{z) 
dz 



V{z) V{-z) 



(20) 



we can rewrite the metric (^ as: 
ds^ = 



N^dsfa) + AN^e^^dzdz 



N'dsf,^ + 



did+1) PiN)^^,^ 



2A P{N) de^ 
d{d+l) iV^ ^ 



2A 

Finally, we can change variables from t to N using (|16|) and write the metric as: 



(21) 
(22) 



ds' 



NV,, + ^dN^ + 2A PiN)dO^ 



2A P{N) 



d{d+l) N^-^ (3^ 



T^^This requirement can be dropped for d = 1. 



(23) 
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where now varies between the two roots, A^i and N2 of P{N). The region of variation of 
the angle 6 depends on the precise form of V{z) and the choice of coefficient (3; (3 in (pOf ) is an 
arbitrary constant that can be reabsorbed in redefinition of 6^ 

Clearly the metric ( ^31) has singularities at = Ni,N2. To find what these singularities 
imply for the physical interpretation of the solution, it is convenient to consider again the metric 
in the form (|2l|) , (p!9|). Clearly the zeros of V{z) are singular points of the metric. Under the 
map z ^ t (^) these zeros are mapped to t = ±00, which, finally, are mapped to the roots of 
the polynomial P{N), Ni^2- 

Now recall that the geodesic equation — the equation of motion of the branes (|^) has also to 
be obeyed (for example, if this equation is not obeyed, an initially static brane will accelerate). 
In the static gauge, = a^,X^^ = 0, (|) implies that y/gT^f^i^g^'^ = at the positions of 
the branes (note that in the static gauge the induced metric is equal to the the metric of the 
embedding space restricted to the brane). For our ansatz (|^) the geodesic equation in the static 
gauge amounts to the condition 

^id-2)x-2K = N'^-^V = , (24) 

which should hold at the positions of the branes. The positions of the branes, therefore, 
correspond to zeros of the holomorphic function V ([TT|), and hence to A^^i 2 in the coordinates 

(H). 

In the following, consider the simple ansatz for V{z), V{z) = z/c; note that since t = 
2A(c log z + c* log z) / {d + l)d, single valuedness of the map t{z, z) requires c to be real. To see 
what kind of singularity this is, consider the /x, v components of the Einstein equations in the 
static gauge: 

sfg (R,u-\g,uR + 9,uI^ = S\y - X,) ~g'/' ~g,, , (25) 

where the hats indicate that the curvatures and metric are those of eq. (|^. The metric in the 
2;-coordinates ( ^I]) can be written in the form: 

rfs^ = dsf^) + e^^ dz dz , (26) 

where 

Near the singularity, equation (p5|) becomes: 

- 1 e'" R2 = -ddK = - - Xa) , (28) 

where the equality holds up to terms that do not contain delta functions. To calculate the 
function k (^7|) near z = 0, 00 note that these two points are mapped by ( pO] ) to t = ±00 
(whether 2; = is mapped to t = +00 or —00 depends on the signs of A and c). These points, 

^^In d = 1, a = 0, P{N) = N'^ + 6, so upon changing variables N — |&|^^^ cosil in ( p3| ) we recover the solution 
obtained in |I| by explicitly solving the equations of motion. 



P{N) 2A 
|z|2 d{d + l) 



(27) 
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on the other hand, correspond to the roots of the polynomial P, A^i and A''2. Near z = 

(assume, for concreteness, that z = is mapped to = A^^i), the dependence of on t is easy 
to find upon solving eqn. (|1^) by keeping the most singular term: 

\N{z)-Ni \ ~ , (29) 



*i = ^^'-^ ■ (30) 



where 

2Ac P'{Ni 
d{d + l) Nt 

Therefore, near z = 0, the function k from the metric (^) behaves as: 

K ^ - logl^l^^^i . (31) 

Upon comparison with eq. (|28|) , noting that ddn = — 27r(l — 6i)6'^{z) we can now find the 
tension of the brane at z = 0, 

f^ = MUn (l-(5i) . (32) 

We can repeat the same analysis near z = oo mapped to the other root, A^ = A^2-Q It is 
convenient to change variables to u = 1/z. Note that the metric ( ^Bf ) (and the function k) has 
the same form in the u variables. The behavior of A^ near u = from eq. (0) is: 

\N{u) - N2\ -\u\~^^\ (33) 



where 

2Ac P'{N2 



d{d + l) Nt^ 



(34) 



"2 

Consequently, k ~ — log and the tension of the brane near z = cxo is then 

= MHtc (1 + 52) . (35) 

We note the conditions that < 5i < 1 and — 1 < ^2 < 0. They follow from requiring 
positivity of the brane tensions (^2 > — 1 and 6i < 1). On the other hand, the distance between 
the two branes should be finite along any path in the z-plane, which requires 6i > and ^2 < 0; 
this follows from considering the interval (|2l|) in the 2;-coordinates and requiring integrability 
of ds near z = and z = oo.Q Finally, positivity of the metric (^) of the transverse space also 



requires that AP{N)N^~'^ be positive as A^ varies over the interval Ni,N2. These restrictions 
impose constraints on the various parameters in the action and the solution; these will be 
considered for our d = 4 example in Section 4. 

We can now deduce a general formula for the Planck scale and cosmological constant on the 
a-th brane. The Einstein term in d-dimensional effective action will be: 

-M'' J d'^x dN J ^-det(Ar2^(,)) R^,^(N^g,) (36) 
' d^x^ R^d){9d) j^' N^-^dN (3-^ j de . 

*^The leftmost equation in ( |2^ ) implies that the roots of P{N) are stationary points of A and the equation of 
motion of the branes is also satisfied at 
'^^Equivalently, > 0, ^2 < ensure that the map \z\ ^ N is nondegenerate near the ends of the interval, as 



follows from the expressions (29, p3| 
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Let the warp factor on the a-th brane be N(^a)- The physical metric there is the induced metric 
^{a)9{d)- Hence the d-dimensional Planck constant — the coefficient in front of the curvature 
term in the action — will be: 




where the 6 integral is over the appropriate range that follows from (^0]) and depends on the 
choice of V and (3. Since the induced metric on the a-th brane is N^^^g(^d), and since g(^d) obeys 
the equation (H), the cosmological constant on the a-th brane, denoted by A(^d),(a), is 

A(d),(a) = iVf,? a . (38) 

Similarly, one can derive an expression for the proper distance between the branes at A^i and 
N2, which we will present for the case d = 4 in the following section. 



4 An explicit example and some speculations. 

As an application of the general analysis of the previous section, let us now consider the 
physically interesting case d = 4. The polynomial P{N) then is of degree 5. It is easy to see 
that if the coefficients a, b in P{N) obey^ 

a > and < 6 < .18 a^/^ (39) 

then there are two positive roots A^^i and N2 of P{N)] in addition, P{N) < as iV varies 
between the two roots, while P'{Ni) < and P'{N2) > 0. Since P{N) is negative in the range 
of variation of N, the metric (^) is positive definite only for negative A, corresponding to de 
Sitter space in our convention.^ Note that since a ~ a/A > 0, this implies that a < as well 
and the four dimensional metric on the branes is also de Sitter (also note that a solution with 
a = does not exist — then P{N) has a single real root and the equation ( |T6|) can not be solved 
for the desired interval of values of t). 

We need to satisfy the conditions Si > for the root Ni which is the image of 2; = and 
6j < for the other root Nj (the image of z = 00), so that the proper distance between the 
two roots is finite. Since 6i ~ AcP'{Ni)/Nf, (^, Q), for our choice of two positive roots, 
< A^i < N2, we find that sign di = sign c and sign 62 = —sign c (since A < 0). On the other 
hand, since t ~ Ac log \z\, if c > 0, then z = is mapped to t = +00 (and, since f{N) < 0, the 
smaller root A^^i), while z = cxd is mapped to t = —00 and, hence, the larger root A^^i. Then the 
finite proper distance requirement is satisfied by the choice c > 0, yielding 61 > 0, 62 < 0. The 
coordinate 6 ( pOf ) changes then from (0, 2ttc) (we chose (3 = —1/2). Expressing the parameters 



^^More precisely, the numerical constant on the r.h.s. is (3/5)^/^ — (3/5)^/^. 

"^^^This conclusion is always true: in the case of two negative roots of P{N), which can be achieved by adjusting 
the value of 6, P > between the roots, but < and positivity of the metric also requires negative A. 
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5i^2 through the tensions of the branes we obtain: 



lAlc \5Nf-3aNf\ ff 



10 Nf AttM^ 

\A\c 5N^-3aNl ^ /| 



The induced metric on the i-th brane has 



A-TTf 

A4, = Nr'a and M|,,, = — Nr^ {N^ - ATf) , (41) 
while the proper distance between the branes is: 

\ A Jn, \N^-am + b\y^ ■ ^ ' 

There are four parameters in the action with which we started: the 6 dimensional gravita- 
tional constant, M, the cosmological constant, A, and the tensions of the two branes, /i,/2- 
The solution involves three additional dimensionful parameters, the parameter a of the ansatz 
for the d-dimensional metric (H), the integration constant e (0), and c — the "strength" of the 
zero of V{z). The equations (P0[ ) expressing the tensions for the branes through the other pa- 
rameters of the solutions represent two constraints on the seven parameters M, A, /i, /2, e, a, c. 
We can use them to eliminate the tensions /i^2 from the list of our parameters. The only im- 
portance of eqns. (|40|) then is that their left hand sides have to be smaller than 1, to ensure 
positivity of the brane tensions. 

From the four dimensional effective theory point of view the quantities of interest are the 
four dimensional Planck constant, Mpi and the cosmological constant, A4. Of interest are 
also the mass of Kaluza-Klein excitations as well as the distance scale where a gravitational 
experiment on the brane will reveal deviations from Newton's law. The last two quantities 
depend, of course on the size of the extra dimensions (and hence on the distance R (|i2|)). 



However, due to the warped geometry, there can be also nontrivial dependence on the warp 
factors; we leave a detailed investigation of this issue for future work. 

To get an idea as to what the effect of the nontrivial warp factors might be, consider first 
the case where the parameters a, b in P{N) (p!7|), (|TBp satisfy < 6 << a^/^. The advantage of 
this limit is that the roots can be easily approximated by: 

Note that the case we are considering corresponds to N2 ^ A^^i — the warp factors on the two 
branes can be (upon adjusting the parameters, of course) very different. 



The conditions of positive tensions (^0]) then imply 



and \A\ c a'/' < 1 . (44) 



10 \b 
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In the case we are discussing, the proper distance (^21) between the branes is easily seen to be 
approximately 

Rr^^ . (45) 



A 

The induced metrics on the two branes at A^^i and are de Sitter with cosmological constants: 



|A(4),i| = \a\ 



\a\ I — 



j and |A(4),: 



|A| 



while the corresponding Planck constants are: 





a 


13/6 


|A|3/2 


e 


2/3 



and M 



PL2 



M^c a"^Nz^ ~ M^c 



1/2 



|A|V2 



(46) 



(47) 



The vacuum energy density on the brane is of order Mp; A(4). If one of the branes described 
the physical universe, the parameters should be tuned so that A(4)^j/M|,; ^ ~ 10^^^° (or to a value 
smaller than that; note also that this ratio is the same on both branes, since the warp factors 
cancel between numerator and denominator). It is not our objective to solve the cosmological 
constant problem; we will only ask whether our parameter space allows for at least a fine-tuned 
situation, subject to the constraints imposed by the existence of a solution. 

Tuning the cosmological constant to the observed small value, therefore, requires that 



A2" 



(|A|c) ^ 10 



120 



(4J 



It is convenient to work in terms of the dimensionless ratios a, 6, (M^|A|~^), (|A|c), and keep 
a as the only dimensionful parameter. 

Consider first the case where the observed universe is at the first brane. Then, 



M; 



PLl 



(|A|c) 



\OL\ 



b) 



2/3 



Hence, 



Then 



2/3 



10 



120 



a\2/3 



10^^ eV 



a J 



(49) 



1 

R 



|A| 



1/2 



~ \a 



10~33 eV 



This is a very large distance, unless b ^ a^/^, which violates (^). If the observed universe is 
on the other brane, we obtain instead |a;| ~ alO^^^el^^ and 1/i? ~ 10~'^'^e\^ — also a large value 
for R. 

Therefore, if the deviations from Newton's law reveal themselves to an experimentalist on 
the brane at a distance scale of order R, then we can not accommodate the observed universe 
on one of our branes. If, on the other hand, the relevant scale of deviations is NiR, then it 
is possible to fine-tune the parameters to obtain a value for N^R of order a millimeter, and 
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hence accommodate the observed world in our setup (admittedly with enormous fine-tuning of 
parameters!). 

The above considerations apply also in a more general situation, where eq. ( ^S]) is replaced 

by: 

^(lAle)M^.lO-. (50) 

while (1^ ) becomes: 

|a| ~ N'l 10-^^e1/2 , (51) 

provided that the observed universe is identified with the brane at Ni. But then the inverse 
radius becomes: 

i ^ ^ / 10--eV . (52) 

where / is the integral 

.iV2 N^I'^dN 
In, \N^-aN^ + h\y^ ' ^^^^ 

which is trivial in the limit 6 — > considered before. 

Assuming that in our geometrical setup the deviations of gravity from Newton's law will 
show up, to an observer on "our" brane, at a distance scale of order i?, an acceptable value of 
1/R > lO^^eV requires Nila^^/'^ > 10'^'^ (which appears hard to achieve, similar to the example 
given above). If the deviations show up at a scale NiR, then the requirement is la'^/"^ > 10^° 
(which is achievable via fine tuning). We leave the investigation of this issue for future work. 
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